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The transport coefficients for a gas of smooth, inelastic hard spheres are 

obtained from the Boltzmann equation in the form of Green-Kubo relations. 
The associated time correlation functions are not simply those constructed 
from the fluxes of conserved densities. Instead, fluxes constructed from the 
reference local homogeneous distribution occur as well. The analysis exposes 
potential problems associated with a straightforward extension of linear re- 
sponse methods to granular systems. 

I. INTRODUCTION 

The derivation of macroscopic equations (e.g., hydrodynamics) and associated transport 
coefficients from kinetic theory is limited in practice to weakly coupled systems (low den- 
sity gases, ideal plasmas, anharmonic crystals). The application of formal methods from 
nonequilibrium statistical mechanics to this problem forty years ago provided the generality 
missing in kinetic theory These methods lead to the expected macroscopic dynamics 
with formally exact expressions for the transport coefficients known as Green-Kubo (GK) 
expressions. The latter are time integrals of equilibrium time correlation functions for the 
fluxes associated with conserved densities. The simplest example is the diffusion coefficient 
D for an impurity in a simple fluid 
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D= lim ~/V(v(*)-v), (1) 
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where d is the dimension of the system, v is the impurity velocity, and the brackets denote 
an equilibrium ensemble average. Similar expressions apply for the viscosities, in terms of 
the autocorrelation function for the momentum flux, and the thermal conductivity, in terms 
of the autocorrelation function for the energy flux. The GK expressions provide an exact 
starting point for analysis and modeling of strongly coupled systems. In particular, they 
have stimulated extensive studies of transport via molecular dynamics simulation of these 
equilibrium time correlation functions. Much of the present knowledge about transport 
in strongly coupled systems (dense gases, liquids, plasmas, solids) derives from analysis of 
appropriate GK expressions. 

Currently, there is great interest in the foundations of a fluid dynamics for granular 
media. Kinetic theory and statistical mechanical methods have been applied to idealized 
granular fluids comprised of smooth hard spheres interacting with inelastic collisions |2| . 
The derivation of the corresponding Navier-Stokes level equations from kinetic theory at 
low and moderate densities has been given recently, with transport coefficients expressed as 
functions of the restitution coefficient . The results are in good agreement with selected 
tests using molecular dynamics and Monte Carlo simulations ||]. It is tempting to expect 
that some form of the more general GK expressions should apply for granular media as 
well fl]T0|j. The primary difference from normal fluids is the absence of a reference stable, 



stationary Gibbs state in terms of which the time correlation functions would be defined. 
Instead, the corresponding homogeneous state for an undriven system has a time dependence 
due to loss of energy on collisions ( "cooling" ) . This is referred to as the homogeneous cooling 
state (HCS) [|11|,|1^]. Recently, analysis of impurity diffusion in a granular fluid has led to a 
GK expression similar to ([!]) with the velocity autocorrelation function defined in the HCS 
]. A related analysis of impurity mobility also gives a GK expression for the mobility 



coefficient |14|], although not simply in terms of the velocity autocorrelation function as is 



the case for normal fluids. The difference is due to the replacement of the Gibbs state by 
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the HCS. This is an indication that translation of linear response methods for normal fluids 
to granular fluids requires some care. 

The objective here is to derive GK expressions for a granular gas based on the Boltzmann 
equation. This may seem redundant since the transport coefficients are already known by 
other methods 0. However, the analysis exposes the form of the GK expressions that should 
be expected from more general nonequilibrium statistical mechanical methods. As with the 
mobility, it is found that the time correlation functions for the transport coefficients are 
not simply those for the conserved densities. Instead, they are generated from the HCS 
distribution in a manner required by internal consistency for the solution to the Boltzmann 
equation. The analysis also shows that the HCS correlation functions can be put in the form 
of stationary state averages using appropriate dimensionless variables, including collision 
number instead of real time. 

It is a pleasure to dedicate this work to Bob Dorfman, colleague and mentor, who has 
taught us all so much about hard sphere transport (among many other subjects). 

II. BOLTZMANN EQUATION AND HYDRODYNAMICS 

The simplest model for a granular fluid is a system of N smooth hard spheres or disks at 
low density interacting via inelastic collisions. The collisions are characterized by a normal 
restitution coefficient a < 1, where a = 1 is the elastic limit. An accurate description of 
this gas can be obtained from a formal density expansion of the BBGKY hierarchy, just 
as for normal gases To leading order in the density, the single particle distribution 

function /(r, v, t) obeys the Boltzmann equation, with collisions modified to account for the 
inelasticity, 

($ + v.V)/ = J[fJ]. (2) 

The detailed form of the inelastic collision operator J is not required for the discussion here, 
beyond the properties required for the macroscopic balance equations: 
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Here p = nT is the low density pressure, n being the number density and T the granular 
temperature (with Boltzman's constant set equal to unity). The functional ([f] is the 
"cooling rate", as will become apparent below. The number density and temperature, as 
well as the flow velocity u are defined in the usual way 
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with V = v — u, the velocity relative to the local flow. The balance equations for these 
fields follow directly from moments of the Boltzmann equation using Eq. (0), 



Dm + nV • U = 0, 



(5) 



D t Ui + {mn)-*~VjPij = 
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D t T + — (PijVjUi + V • q) + T( = 0, 
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where D t = d t + u • V is the material derivative. The pressure tensor and heat flux q 
are linear functionals of / given by 



P ij [r,t\f]=p(v,t)5 ij + J dvD y (V)/(r,v,*), 



(8) 



q[r,t|/] = |dvS(V)/(r,v,t). 



(9) 



Finally, the functions Z)y(V) and S(V) characterizing the irreversible parts of the fluxes are 



A,(V) ee mm - \ V %i), S(V) = (|^ 2 - V. 
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Clearly, the balance equations are not closed until the fluxes are expressed as functional 
of the hydrodynamic fields n, u, and T. This can be accomplished if a solution to the 
Boltzmann equation can be obtained in "normal form" , 

f(r,v,t) = f(v\{yp}), (11) 

{yp{r, t)} <— > Mr, t), u(r, t), T(r, *)} . (12) 

The notation f {v\{yp}) means that the distribution is a functional of the hydrodynamic 
fields yp(r, t), and that its space and time dependence occurs only through that of these fields. 
If such a solution can be found, its use in Eqs. (§) and (^) gives the fluxes as functional 
of the fields. Such fluxes are then referred to as constitutive equations. The constitutive 
equations together with the exact macroscopic balance equations become a closed set of 
hydrodynamic equations. The origin of a hydrodynamic description in this context therefore 
is traced directly to the existence of a normal solution to the Boltzmann equation on some 
length and time scales. In practice, the functional form of this solution is constructed in 
some well-defined approximation, leading to corresponding approximate constitutive and 
hydrodynamic equations. The approximation of interest here is small spatial variations of 
the hydrodynamic fields over distances of the order of the mean free path. For gases with 
elastic collisions, the Chapman-Enskog method generates the normal solution perturbatively, 
and to leading order the constitutive equations are Newton's viscosity law and Fourier's heat 



law |lq| . The hydrodynamic equations become the Navier-Stokes equations. Application of 
this method to granular gases leads to similar results J7j. The analysis is modified here to 
obtain the same results in an equivalent form with the transport coefficients represented by 
GK expressions. 



III. SMALL GRADIENT SOLUTION 



The first step in constructing the desired solution to the Boltzmann equation is to express 
it in terms of a reference distribution f(°\ 
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/(r,v,t) = /(°)(v|{ y/3 }) + / (1, (r,v,t). 



(13) 



The reference distribution is restricted to be normal and to give the exact moments of Eq. 
(1), i.e. 
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The Boltzmann equation then becomes 



with the definition 
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(dt + v • V+L) /« - /CD] = J[/(°), /(o)] _ ( A + V ■ V) / 
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The left side of Eq. ( |i~5|) generates the dynamics of The solution of interest is such 
that /W should be proportional to the gradients of the hydrodynamic fields, since 
provides their local values through the constraint (0). This requires that the right side 
of flT5|) be proportional to the gradients of the fields. Evaluation of (D t + V • V) using 
the macroscopic balance equations gives terms proportional to the gradients, except for 
the contribution from the cooling rate in (0). In addition, the nonvanishing contribution 
J[f(°\ f^] also is not of first order in the gradients. Therefore, the reference distribution is 
finally characterized by the condition that these non-gradient terms should vanish 



J[/ (0) ,/ (0) ]+TC (0) ^ = 0, 
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This is an essential point in the analysis, and it is the origin of differences from the case of 
elastic collisions, as is discussed below. The distribution function is the analogue of the 
local Maxwellian for the case of elastic collisions, but the solution to Eq. ( |TTD for a < 1 is 
not the Maxwellian. Symmetry and dimensional analysis requires that f(°' have the scaling 
form 



/ (0) (v|{^}) = V(i)/ (0)> (V*) , V* = V/v (t), (18) 

where Vo(t) = [2T (t) /m] 1 ^ 2 is the "thermal" velocity and is an isotropic function of V*. 
Then Eq. (JT7|) is equivalent to 

J[/ (0) ,/ (0) ]-ic (0) ^-(v/(°))=0. (19) 

Consider the initial condition /^(0) = 0. This is a physically interesting case for hydro- 
dynamics since /^(0) does not contribute to the initial value of the hydrodynamic fields, 
due to (|14D, and /^(0) is entirely determined by the hydrodynamic initial values. In this 
case, is manifestly proportional to the gradients of the fields. Retaining only terms of 
first order in these gradients, the Boltzmann equation becomes 

(d t + L)fW = -fW* p -Vy p , (20) 

with the definitions 

*^( V +^ (0) )' < 21 > 

^-(^V.^-^)ln /( o), (22) 

$ 3 = m -i JL i n ;(0) _ X. ( d + V ■ A In / (0) ^ . (23) 

A term from to first order in the gradients does not occur since it vanishes. This follows 
from the fact that ( is a scalar and can only be proportional to V • u at this order. Since 
$2,y(V) is traceless, there is no such contribution. The solution to Eq. (^) can be written 
in the form 

/(i) = -Fp-Vyp. (24) 

When Eq. ( |24T) is substituted into Eq. (^), the terms from dtWyp are higher order in the 
gradients, except for f3 = 3 which gives a contribution proportional to the cooling rate. The 
coefficients of the gradients are then found to obey the equations 
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IV. SCALING 



The apparent simplicity of Eqs. Q2"5D-(PS|) is misleading since 



L, C*- , and are functions 



of time through their dependence on the hydrodynamic fields. This dependence occurs 
as well for a gas with elastic collisions, but it can neglected in that case when solving 
the equations since it is proportional to higher order gradients. For elastic collisions, the 
temperature gives a time dependence that is of zeroth order in the gradients and cannot be 
neglected. However, it can be removed by a change of variables to dimensionless forms. A 
characteristic length scale is given by the mean free path i and a dimensionless time s is 
defined accordingly by 



For simplicity we will take here I = (ner^ -1 ) -1 , where a is the diameter of the particles, 
omitting a factor that depends on the dimension of the system. Integrating over an interval 
t shows that s is an average number of collisions during that time. Moreover, introduce 
dimensionless functions &*p and F£ as 




(29) 




mv$ +2 (t) 



2n£ 



F*(V*,s). (30) 
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with the definitions 



(d. + r) (Ft - F*) = (*; - *•) , (31) 



^ + ^ + ^C)^,=/ (0) *^,, (32) 



a + £*-±C*)F5 = /<°>*#5. (33) 



£ *= r + ^k V *' (34) 



-A C = \r- (35) 



V (t) ' Uo(*) 

The right sides of Eqs. (j31~l)- (|32l) are independent of s, as are C* and £*. The integration is 
now trivial 

F* (V*, s) = Ft (V*, s) + f ds'e- s ' r f {0) * (*J - *S) , (36) 

JO 

(V*, s) = f d S ' e - s '( r+ ^*)/(°)*^, (37) 
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F* 



(V*,s)= I* ds'e- s '( T -3 c )f( 0) *&t (38) 

JO 



In the derivation of the above expressions, the assumed initial condition, implying that all 
the F^ vanish for s = 0, has been used. 



V. CONSTITUTIVE EQUATIONS 

The pressure tensor and heat flux are determined from their definitions in Eqs. (§) and 
. Substitution of Eqs. flT3]) and ( PH ) leads to results valid to first order in the gradients, 

Pij [r, t\f] = p(r, t)<% - | dvDyWFp-Vyp, (39) 
q[r,t|/] = - JdvS(V)F p -Vy p . (40) 
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Using fluid symmetry, these expressions reduce to 



(41) 



q = -kVT - fiVn. 



(42) 



These are the expected Navier-Stokes constitutive equations, except for the additional term 
in the heat flux proportional to Vn. It will be seen below that this new term is due entirely 
to the deviation of from the Maxwellian. The expressions for the transport coefficients 
are identified as 
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(43) 



(44) 



(45) 



The dimensionless forms of S and Da a are 



s* = ( v* 2 - d -^- 



D*j = V*V* - ^V* 2 6 i:j . 



To put these expressions in the desired Green-Kubo form, define the adjoint operators 
L* and C* by 

( L'^ 







j d\*X(V*) 







Y(V*) = - J dV*Y(V*) 



XV* 
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(46) 



for arbitrary functions X and Y . From Eq. ( p4|) it is trivially seen that 

1 d 

C* = L* + -C*V* . 

2 s dV* 



(47) 



The expression of L* can be obtained by using the properties of the Boltzmann collision 
operator, but it will be not needed here. 
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Also, define "correlation functions" by 

(XY) = J d\*f^*(V*)X{V*)Y (V*). (48) 

Then, with these definitions and Eqs. (|36|)-(|38"|), the transport coefficients given by Eqs. 
(|3|)-(|45|) can be rewritten as 

2nm£vo(t) 
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d 2 + d-2 
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[ S ds' {D^s')^) e-V? , (49) 
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d Jo 
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n 



+ ™(&'(SV)-($l-^)>. (51) 

The time dependence of the correlation functions is given by 

X{s) = e sC *X(V*). (52) 

Equations (|5*l"D-(f4*9"D are the Green-Kubo expressions for the transport coefficients of a low 
density granular gas. The spectrum of C* is such that the correlation functions decay to zero 
for s ^> 1. This sets the time scale for hydrodynamics. In this limit, the above expressions 
agree with those obtained by the Chapman-Enskog method @. 

VI. DISCUSSION 

To discuss the GK expressions for a granular gas it is instructive to write the corre- 
sponding results for elastic collisions (a = 1). In that case, /(°) is the local Maxwellian 
and 

*!^0, -ir. <&*^S*(V). (53) 

Then \i — > fi Q = 0, and 



J* ds' (plj(s')D*^ , ^^^[^'(SVl'ST (54) 



d? + d-2Jo \ l3V ' lJ /' d Jo 
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Since the temperature is constant to leading order in the gradients for a = 1, the variable 
s = (vo/£)t is simply proportional to time. These are the standard forms for the GK 
expressions in terms of the autocorrelation functions composed from the "microscopic" fluxes 
S* and D*j. There are several differences that occur for granular gases: 

• The correlation functions are defined as averages over f(°\ This local HCS distribution 
is determined from fll9T) and differs from the Maxwellian for all a < 1. 



• The correlation functions are not constructed from S* and D*j alone. Each is paired 
with another function from the set 

• The time integration is replaced by an integration over the average collision number, s. 
The correlation functions have approximate exponential decay in the variable s rather 
than t. 

• The integrals over s are not controlled solely by the correlation functions. In addition 
there are time dependent factors arising from the change of the temperature over the 
duration of the integral. 

The most surprising among these differences is the replacement of one of the "micro- 
scopic" fluxes S* or D*- by one of the new variables <E»^ determined from through Eqs. 
(pl|)- (|23|) . If /(°) were replaced by the local Maxwellian f^* in the above analysis, the corre- 



lation functions (S*(s') • S*) and \JD*-{s')D*-j would appear in Eqs. (HD-flO). However, this 
would not be a consistent solution to the Boltzmann equation and the transport coefficients 
would not have the correct dependence on a (e.g. fi would vanish for all a). The role of 
f(o) [ s m ore than just a local reference state with the exact moments for n, T, and u. In 
addition, it must be an approximate solution for the dynamics. In dimensionless form, /(°)* 
is a universal function of V* determined from the equation 

J*[f {0> , / (0) *1 - ■ (v7 (0) *) =o, c* = C*[/ (0) 1, (55) 
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without reference to any particular hydrodynamic state. The above analysis shows that 
(|55l) is necessary for a consistent ordering of the solution in terms of gradients. The local 
reference distribution for a given hydrodynamic state follows from this solution according to 



/<°>=n(r,t) 



m 



d/2 



f 
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vv-u(r,t)) 



(56) 



Its dependence on the hydrodynamic variables is determined by the solution to (|55|). 

An estimate for the can be obtained from an approximation to f^* obtained by a 
polynomial expansion CTTT2J, which to leading order is 
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(57) 



(58) 



9 + 24d + (8(2 - 41)a + 30a 2 (l - a) 
This is known to be accurate to within a few percent for all a and V* < 1. In this approxi 



mation the variables become 



2,y 
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(61) 



Linear response methods typically construct the response function for spatial perturba- 
tion of the homogeneous state of a specific type: those that couple only to the microscopic 
conserved densities. This follows from consideration of an initial local equilibrium state 
or maximum entropy state. Then manipulation of the response function for the conserved 
densities using the microscopic conservation laws leads to GK expressions in terms of the 
microscopic fluxes. However, it is clear that expansion of (|56|) for small spatial perturbations 
does not give a representation in terms of the conserved densities. Thus, a straightforward 
application of linear response can lead to formal results with internal inconsistencies. One 
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such inconsistency can be a failure of the solubility conditions for the linear integral equa- 



tions associated with (|31| ) - fl32|). In the GK form these conditions translate to requirements 
that the fluxes be orthogonal to invariants of the dynamics. These issues and a derivation 
of the GK expressions from the Liouville equation will be given elsewhere. 
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